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STRUCTURAL STABILITY AUGMENTATION SYSTEM DESIGN USING BODEDIRECT: 

A QUICK AND ACCURATE APPROACH 

By 

T. J. Goslln and J. K. Ho 
Boeing Commercial Airplane Company 
Seattle, Washington 


ABSTRACT 


Both aft and forward sections of long slender aircraft suffer from lateral accelerations 
in turbulence. Some of these accelerations can be attributed to flexible body bending of 
the airplane in the 1-6 Hz region. Given an accurate flexible body model of the aircraft, 
a control law can be employed using lateral acceleration feedback loop augmented with 
the basic dutch roll yaw damper system to actively damp out these modes. 

In this paper a methodology will be presented for a modal suppression control law 
design using flight test data Instead of mathematical models to obtain the required gain 
and phase information about the flexible airplane. This approach will be referred to as 
BODEDIRECT. The purpose of the BODEDIRECT program is to provide a method of 
analyzing the modal phase relationships measured directly from the airplane. These 
measurements can be achieved with a frequency sweep at the control surface input 
while measuring the outputs of interest. The measured "Bode-models" can be used 
directly for analysis in the frequency domain, and for control law design. Besides 
providing a more accurate representation for the system Inputs and outputs of interest, 
this method is quick and relatively Inexpensive. 

To dale, the BODEDIRECT program has been tested and verified for computational 
Integrity. Its capabilities include calculation of series, parallel and loop closure 
connections between Bode-model representations. System PSD, together with gain and 
phase maiglns of stability may be calculated for successive loop closures of multl- 
inpul/multi-output systems. Current plans Include extensive flight testing to obtain a 
Bode model representation of a commercial aircraft for design of a structural stability 
augmentation system. 

In addition to the BODEDIRECT approach, an Indirect approach using flight test data to 
derive a mathematical mode for analysis using a Transfer Function Matching Routine 
will be presented along with the strengths and weaknesses of each approach. 
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BOEING COMMERCIAL AIRPLANES 


THE RIDE QUALITY PROBLEM 
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THE MODES WHICH MUST BE ACTIVELY CONTROLLED INCLUDE 


BASIC YAW DAMPER : 

The Dutch Roll Mode 

Spiral Mode for turn coordination 
MODAL SUPPRESSION SYSTEM ; 

1st and 2nd Body Bending Modes 

All other modes must be passively controlled or have 
sufficient gain and phase margins. 


Note : This discussion is limited to the sysnthesis of a modal 
suppression system using the BODEDIRECT program. It is assumed 
that a good basic yaw damper already exists. 
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TRADITIONAL METHOD 


THEORETICAL 

ANALYSIS 


WIND TUNNEL 
AERODYNAMIC 
DATA 


FLIGHT TEST 
DATA 


LARGE STATE - SPACE 
STRUCTURAL MODELS 


NUMBER \J BODE PLOTS 
CRUNCHING X| AND PSD 


CONTROL LAW 
DESIGN 


FIGURE 1 A 


BODEDIRECT 



FIGURE IB 
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Mathematical Models BODEDIRECT 
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1) Require large computing D Qulck and inexpensive (work 

budgets (main frame computer) station environment) 



BODE - MODEL GENERATION 


c 5 

« cd 
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Bode-models are measured for all inputs and 
outputs of interest for the free-airplane 




MODAL SUPPRESSION SYSTEM 
CONTROL LAW DESIGN 



For control law design, we determine the 
modal phase relationship in the feedback 
loop such that a compensator can be synthesized 
to bring the feedback signal into phase with 
the desired mode(s). 
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Frequency points may be loaded or generated in log or linear scale. 



Similarly bode-models may either be loaded from a file 
in terms of magnitude and phase or generated from a 
transfer function. 

Example: 

‘LOAD 'modelname' 

‘GENERATE 'modelname’ 

transfer 
function \ 
coefficients! 

In addition, a bode-model may be extrapolated from a previous 
bode-model or a set of lab test data using the command: 


‘CREATE BODE-MODEL 'modelname' 

data I 

points to < 
extrapolate 1 


This allows quick bode-model generation from a few test data points. 
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B0DE-M0DEL1 


B0DE-M0DEL2 




1 




GAIN 

► 

I PHASE 



PHASE 

u 

r 


‘CONNECT IN SERIES : 

MAG = GAIN * MAGI * MAG2 
PHASE - PHASE1 + PHASE2 
IF( GAIN .LT. 0 ) PHASE - PHASE + 180° 

‘CONNECT IN PARALELL 

MAG-V<COM 2 + YCOM 2 
PHASE - TAN - 1 (YCOM/XCOM) 

where 

XCOM - GAIN1 * MAGI * COS (PHASE1) 

+ GAIN2 * MAG2 * COS (PHASE2) 
YCOM - GAIN1 * MAGI * SIN (PHASE1) 

+ GAIN2 * MAG2 * SIN (PHASE2) 
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"RESOLVE FEEDBACK JUNCTION 


OUT 


MAG * l/V XCOM 2 + YCOM 2 
FBAC 

PHASE - -TAN ' }YCOM/XCOM) 

FBAC 

where 

XCOM - 1.0 - MAG (G*H) * COS (PHAS (G*H)) * GAIN 
YCOM - -MAG (G*H) * SIN (PHAS (G*H))* GAIN 



"RESOLVE CLOSED-LOOP SYSTEM 
CLOSED-LOOP - FBAC * G (SERIES CONNECTION) 

"COMPUTE BODE-LOCUS 

(This routine computes a locus of bode-plots for a range 
of specified closed-loop gains. The user provides G,H, and 
the range of gain values.) 


"COMPUTE GAIN AND PHASE MARGINS 

(The gain and phase margins are determined from the broken loop 
system bode plot. For good stability the gain should be -6 db or 
less at the ±180° crossings. Additionally, the phase should be ±45 
away from ±180° at the 0 db crossing. For clarity, all phase 
values are scaled between 0 and -360 in all bode-model plots. 
Uncertainty values for gain(db) and phase(deg) may be declared by 
the user. The program will warn of near 0 db or 180 deg crossings 
based on these uncertainties.) 


I 
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EXAMPLE USING SIMPLE (YAW RATE ONLY ) BASIC YAW DAMPER 
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•RESOLVE OPEN-LOOP SYSTEM 


THIS ROUTINE WILL RESOLVE THE OPEN LOOP BODE-MODELS 
(G1 , G2, G3...) FOR A MULTILOOP SYSTEM IF THE CLOSED-LOOP AND 
FEEDBACK BODE-MODELS ARE KNOWN. A GAUSS-JORDON 
SOLUTION IS REPEATED FOR EACH FREQUENCY POINT TO SOLVE 
THE SET OF SIMULTANEOUS EQUATIONS. 


U 



THIS ROUTINE MAY BE USEFUL WHEN THE FLIGHT TEST DATA IS 
OBTAINED WITH THE BASIC YAW DAMPER ON AND THE EFFECTS OF 
THE YAW DAMPER ON STRUCTURAL MODES ARE THE DESIRED 
INFORMATION. 
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BODEDIRECT INCLUDES NONLINEAR CAPABILITIES 


which incorporates a describing function (DF) analysis 
for the following : 


a) hysteresis 

b) rate or position limits 

c) deadzones 

d) variable gain 


common nonlinearities 


in control actuation 


(rudder) systems 


Using this analysis, different bode-models may be 

generated as a function of input amplitude. Using 

polar plots for the describing function a prediction 

can be made if limit cycling will occur and at what 

frequency and amplitude. Notice that BODEDIRECT 

may use actual lab test data for the rudder system 

if desired over the theoretical model. 
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SPECTRUM 

(G 2 ) 


CRUISE CONDITION 
SQUARE WAVE SWEEP DATA 
LATERAL ACCELERATION AT PILOT STATION 
BODEDIRECT PERFORMANCE PREDICTION 
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<940 


CRUISE CONDITION 
SQUARE WAVE SWEEP DATA 
LATERAL ACCELERATION AT AFT STATION 
BODEDIRECT PERFORMANCE PREDICTION 



POWER 

SPECTRUM 

CG 2 ) 


CRUISE CONDITION 
SINE WAVE SWEEP DATA 
LATERAL ACCELERATION AT AFT STATION 
BODEDIRECT PERFORMANCE PREDICTION 
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MSS OFF 
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PSD 


DESCENT CONDITION 
SQUARE WAVE SWEEP DATA 
LATERAL ACCELERATION AT PILOT STATION 
BODEDIRECT PERFORMANCE PREDICTION 
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PSD 


DESCENT CONDITION 
SQUARE WAVE SWEEP DATA 
LATERAL ACCELERATION AT AFT STATION 
BODEDIRECT PERFORMANCE PREDICTION 
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DESCENT CONDITION 
SINE WAVE SWEEP DATA 
LATERAL ACCELERATION AT PILOT STA1 
BODEDIRECT PERFORMANCE PREDICT!! 
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PSD 


DESCENT CONDITION 
SINE WAVE SWEEP DATA 
LATERAL ACCELERATION AT AFT STATION 
BODEDIRECT PERFORMANCE PREDICTION 
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